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We have observed that the tunneling magnetoconductance between two-dimensional (2D) electron 
gases formed at nominally identical InAs-AlSb interfaces most often exhibits two sets of Shubnikov- 
de Haas oscillations with almost the same frequency. This result is explained quantitatively with 
a model of the conductance in which the 2D gases have different densities and can tunnel between 
Landau levels with different quantum indices. When the epitaxial growth conditions of the interfaces 
are optimized, the zero-bias magnetoconductance shows a single set of oscillations, thus proving that 
the asymmetry between the two electron gases can be eliminated. 
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The unusual energy-band alignment of heterostruc- 
tures that combine InAs and GaSb (or their alloys) 
has made them attractive for electronic devices like 
interband-tunncling diodes and quantum-cascade lasers. 
Lately, those structures have also emerged as one of the 
leading candidates for "spin" devices, in which the elec- 
tron's spin states rather than its charge determine the 
operation of the device, be it a transistor^ or a resonant- 
tunneling diode 

In a III-V compound-semiconductor heterostructure 
with an asymmetric potential profile, the two spin sub- 
bands of the conduction band are split for non-zero in- 
plane wave vectors in the absence of a magnetic field. 
This splitting, frequently called Rashba splitting)^ is 
particularly large (several meV) in InAs because its band- 
structure parameters favor a large spin-orbit-coupling co- 
efficient, and may lead to two populations of electrons 
with opposite spin orientations. 

On the long road to practical spin devices, one of the 
first steps has been to establish unequivocally the pres- 
ence of Rashba splitting in a given heterostructure, to 
determine its amount, and, if possible, to control it. 
The effect has manifested itself as a beating pattern 
in the Shubnikov-de Haas (SdH) oscillations of the in- 
plane magnetoresistancei^ Sometimes, an external elec- 
tric field perpendicular to the layers has been used to 
further modify the relative carrier densities^ However, 
some experimental results on the Rashba splitting have 
been inconclusive. For instance, Brosig et al. did not 
observe any beating in the SdH oscillations of asymmet- 
ric InAs-AlSb quantum wells with various carrier den- 
sities^ and Heida et al. found no significant effect on 
the beating pattern by a gate voltage^ Moreover, the 
presence of a beating pattern in magnetoresistance os- 
cillations may signal a phenomenon quite different from 
spin splitting, for example mixing of the first-subband 
series with magneto-intcrband oscillations^ 

Since some of the proposed spin devices rely on electron 
tunneling across a heterostructure, it is equally impor- 



tant to sort out effects in vertical transport that might be 
construed as evidence of spin splittingiii Recently we ob- 
served a beating pattern in the zero-bias tunneling mag- 
netoconductance of an InAs-AlSb-GaSb-AlSb-InAs het- 
erostructure with a potential profile adequate for a spin- 
filter device^ However, that pattern was explained not 
as a consequence of Rashba splitting but in terms of tun- 
neling between two-dimensional electron gases (2DEGs) 
with slightly different densities^ This conclusion was 
based on an analysis that used a very simple magneto- 
tunneling model, in which only one Landau level at a 
time participated in the tunneling process. 

In this paper, we report a systematic magnetotunnel- 
ing study using a variety of samples with different thick- 
ness for the GaSb well and the AlSb barriers, and employ- 
ing a more elaborate transport model, which allows for 
tunneling involving several (energy-broadened) Landau 
levels simultaneously and in which the Landau-level in- 
dex may not be conserved in the process. This study not 
only confirms the generality of our initial conclusion but 
also shows how the epitaxial growth of the heterostruc- 
tures affects dramatically the beating pattern, which dis- 
appears when the materials interfaces are optimized. 

The band profile common to the InAs-AlSb-GaSb- 
AlSb-InAs heterostructures discussed here is depicted 
in the central inset of Fig. 1, and in greater detail in 
Ref. 0. Since the top of the valence band of GaSb is 
higher in energy than the bottom of the InAs conduc- 
tion band, two-dimensional (2D) electrons accumulate in 
quasi-triangular wells formed at each of the two InAs 
interfaces. Ideally, the number of 2D electrons in each 
interface should be the same, and the sum of the two 
equal to the number of 2D holes left in the central GaSb 
valence-band well. 

To illustrate the similarities and differences among 
the many structures we have studied, we focus on four 
heterostructures prepared by molecular beam epitaxy 
(MBE) at different institutions and in a time span of 
almost a decade: sample A, with GaSb and AlSb thick- 
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nesses of 75 A and 25 A, respectively; sample B, with 
corresponding 60 A and 34 A widths; and samples C and 
D, both with 82 A and 31 A GaSb and AlSb regions, 
respectively^ Samples A and B were grown a few weeks 
apart but several years earlier than samples C and D, 
which were grown on the same day, and as explained 
later, differed from each other only on the procedure fol- 
lowed to prepare one of the InAs-AlSb interfaces. 
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FIG. 1: Shubnikov - de Haas oscillations of various InAs- 
AlSb-GaSb-AlSb-InAs heterostructures described in the text: 
(a) sample A at T=4.2 K, (b) sample B at T=1.7 K, (c) 
sample C and (d) sample D at T=1.7K. The insert illus- 
trates the potential profile common to all samples, showing 
the InAs conduction band, the GaSb valence band and AlSb 
(not listed) in between. 2D electron (e) and hole (h) gases 
reside in InAs and GaSb layers, respectively. Ef is the Fermi 
energy. 

The central panel, (b), of Fig. ^ shows the zero-bias 
tunneling conductance versus magnetic field (perpendic- 
ular to the layers) for sample B, the focus of our ear- 
lier study>i£ As discussed then, the magnetoconductance 
reveals SdH oscillations with a beating-like pattern be- 
tween 0.5 T and 1.5 T, followed by doublet structures 
(unrelated to spin splitting) with varying relative inten- 
sities up to about 4 T. The conductance of sample A 
[top panel, Fig. 2(a)] exhibits a behavior qualitatively 
similar to that of sample B, although the characteristic 
fields and relative intensities of the doublet oscillations 
are somewhat different. 

In sample C [Fig.^c)] there are also a beating-like pat- 
tern and a doublet region, but they are "squeezed" to a 
narrow field interval, between 0.4 T and 1.2 T, so some of 
the features are extremely sharp and easy to miss. Above 
1.2 T a simple set of oscillations with a single period de- 



velops, and Landau-level spin splitting becomes visible 
at 3.5 T. Finally, the oscillatory behavior of the conduc- 
tance of sample D [Fig. Q[d)] is simple throughout the 
entire field range, without any trace of beats or doublets 
(aside from spin splitting above 2 T). 

A fast-Fouricr-transform analysis of the conductance 
oscillations revealed the main frequency component, but 
failed to yield clear evidence of more than one compo- 
nent, except for sample B. This is not too surprising, in 
view of the poor beating pattern observed, particularly 
for samples A and C. When analyzed with the help of a 
simple model for the tunneling magnetoconductance, the 
beating and doublet structures found in sample B were 
attributed to small differences in the carrier densities of 
the 2D electrons at the InAs interfaces^ In the follow- 
ing, we describe a similar approach to analyze the results 
of samples A to C using a model that keeps the essence 
of the previous one but is more general in the treatment 
of tunneling between Landau levels. 

The basic assumption underlying our model is that the 
oscillations in the conductance observed in Fig. ^ are di- 
rectly related to magnetic states of electrons in the InAs 
regions only. This assignment is based on the temper- 
ature dependence of the amplitude of those oscillations, 
from which an effective mass is extracted that agrees with 
the electronic effective mass in InAs. In our model we ac- 
count for the Landau quantization of 2D electrons in both 
InAs accumulation regions but ignore the Landau quan- 
tization of holes in the GaSb quantum well. The hole 
gas is treated simply as a "window" through which elec- 
trons tunnel from one accumulation layer to the other. 
This simplification is reasonable, considering the small 
cyclotron energy of heavy holes (when compared to that 
of electrons in InAs) and the fact that we have not found 
any features in the experimental magnetoconductance 
that could be attributed to Landau-quantized holes. 

With those premises, we model the conductance as due 
to tunneling between 2D electron gases (in general, with 
different electron density) separated by a potential bar- 
rier, first assuming that in-plane momentum is conserved 
throughout the process. We follow Lyo's treatment, in 
which the tunneling conductance is expressed asi^ 
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Jo is the zero-field tunneling integral, g is the Landau 
level degeneracy per spin, and rii is the Landau index 
for the ith layer. /'(C) is the energy derivative of the 
Fermi-Dirac distribution function at energy (. The term 
<5n 2 .ni (= 1 when ri2 = n\ and otherwise) reflects the 
conservation of in-plane momentum in the presence of a 
perpendicular magnetic field. pi n is the density of states 
(DOS) of the n t h Landau level in the ith layer (i = 1,2). 
To approximate Eq. we use a Gaussian density of 
states 
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FIG. 2: Various contributions to the magnetoconductance 
from tunneling channels whose initial and final Landau levels 
have indices that differ by An from to 3, calculated using 
the parameters for sample B (described in the text). 

with a half-width broadening 



where £i n is the energy of the n t h Landau level in the 
ith layer, 77 is a materials-dependent constant, to c is the 
cyclotron frequency, and H is the magnetic field perpen- 
dicular to the interfaces. 

The chemical potential, n, is determined by keeping 
the total carrier density constant 



with 




and A = eH/nh, 7 = 2T/k B T, z = (/k B T, and /j,' = 
fx/k B T. 

In the zero-temperature limit, Eq. Q gets simplified 

to 

G = 4 ^ - Jq Pi»(m)P2w(m " A ) , (5) 
n 

where A = £2 — e± is the energy difference between the 
two subbands at H = 0, and therefore is proportional 
to the 2D carrier difference. Eq. i s a l so simplified 
accordingly. 

With these T = K equations, we can easily calcu- 
late the zero-bias tunneling conductance between two 2D 
electron gases with different carrier densities. The main 
trace in Fig. [5] (curve An = 0) is the result of such a 
calculation, when the following parameters (correspond- 
ing to sample B) are used: m e = 0.027 mo, A = 4 meV, 
N s = 5.7 x 10 11 cm" 2 , and n = 0.8. As seen in Fig. 
even though A ^ the calculated magnetoconductance 
exhibits only one set of oscillations, contrary to the ex- 
perimental result of Fig. mb). The single set in the cal- 
culation is understandable: in Eq. (JSJ the product of 



two density of states for Landau levels whose centers are 
shifted by an energy A but have the same index (and 
therefore the same rate of change with field) still looks 
just like one single density of states^ 

To explain the existence of more than one set of oscil- 
lations, in addition to different carrier densities in the 2D 
accumulation layers one needs to include the possibility 
of tunneling between Landau levels with different quan- 
tum numbers, in other words, to relax the in-plane mo- 
mentum conservation law. Let us consider emitter and 
collector accumulation layers whose zero-field subband 
energies differ by an amount A as a result of different 
carrier densities. At low magnetic fields, that is, when 
A is larger than or comparable to the level broadening 
r and to the cyclotron energy %ui c , the Landau levels 
in one electrode are significantly misaligned with respect 
to those in the other electrode. Then, for a fixed field 
the dominant contribution to the density of states at the 
Fermi level will correspond to one Landau level for the 
emitter and to another for the collector [see Fig. Eta)]. If 
we now allow for tunneling transitions in which An 7^ 
more tunneling "channels" become open in Eq. l(Tj). 

On the other hand, at high fields, when hoj c and T are 
much larger than A, the misalignment between Landau 
levels with the same index in the two electrodes becomes 
relatively insignificant [see Fig. EJb)]. In other words, 
the contribution of the An 7^ channels to the conduc- 
tance is important at low fields but diminishes rapidly as 
the field increases. The situation is illustrated in Fig. [21 
where, in addition to the An = channel, we show the 
individual contributions to the calculated conductance of 
the An = 1 — 3 channels. 
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FIG. 3: An illustration of tunneling betwen Landau levels in 
which the Landau index may not be conserved. Two 2DEGs 
with subband-energy difference A are separated by a barrier. 
Ef is the Fermi energy, n, n' are the Landau indices of the 
right and the left 2DEGs, respectively. The shadowed re- 
gions correspond to filled states. Electrons of Landau index 
n tunnel into (a) n' = n and n -1, or (b) the same index, n, 
whenever the Fermi energy falls into the DOS of that Landau 
level. 

In order to account quantitatively for the experimental 
dependence of the conductance on magnetic field, it is 
important to keep in mind that the contributions of the 
various An channels to the conductance may not all be 
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the same. Their relative weight will depend on details of 
the scattering mechanisms responsible for the violation 
of the An = rule. In the absence of a theory that 
includes those details, we have determined that weight 
empirically from a fit of the experimental conductance 
to an appropriately modified Eq. l|T]). 

Figure 0] compares the experimental and calculated 
conductance for sample B, using adequate probability ra- 
tios for An = 0, 1 and 2. A good fit is obtained when 
the probability is proportional to (1+ | An |) -1 ^ 2 . As 
seen in the figure, for fields approximately below 2 T the 
An = 2 and An = 1 transitions are successively domi- 
nant. At about 2 T, the An = 1 and An = transitions 
have comparable strength and their corresponding oscil- 
lations similar height, as highlighted by the short line 
over the experimental curve. Above 3 T the conductance 
is almost completely governed by the An = transition. 
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FIG. 4: A calculation (lower curve) for the tunneling conduc- 
tance that includes several An channels reproduces well the 
features seen in the experimental result (upper curve) for sam- 
ple B. The numbers 0, 1, and 2 correspond to the difference 
between the indices of initial and final Landau levels, An. 
A short line above the upper curve around 1.8 T indicates 
the change of the dominant contribution to tunneling, from 
An = 1 to An = 0. The insert shows the field dependence of 
this change when a small bias is applied to the sample. 

The "cross-over" magnetic field (about 1.8 T in Fig.0} 
at which the An = 1 and An = contributions to the 
conductance are comparable is closely related to the dif- 
ference between the subbands energies of the two 2D elec- 
tron gases, or, equivalently, to the carrier asymmetry be- 
tween the two gases. This asymmetry is affected by a 
small bias between the electrodes, which slightly changes 
the potential profile along the tunneling direction and the 
carrier density in the accumulation layers but does not 
alter the near-equilibrium condition (insignificant tunnel- 
ing current) of this study. In the insert of Fig. 0] we show 
the value of the cross-over field for several biases, from -2 
mV to +2 mV. With increasing positive bias (of the top 
electrode relative to the substrate electrode) the cross- 
over field increases, indicating an increasing asymmetry. 



TABLE I: Average carrier density, N ave , difference between 
the carrier densities of the two gases, AAf e , and energy dif- 
ference between their two ground subbands, A (using m e = 
0.027 mo), for samples A to D. N ave was determined from 
a direct analysis of the SdH oscillations, while AN e was ob- 
tained from a detailed comparison between the theoretical 
and experimental magnetoconductance oscillations^- 
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On the other hand, with negative bias the field becomes 
smaller and the asymmetry decreases. 

A fit of the calculated conductance to the experimental 
one yields the individual carrier densities in the 2D gases. 
The difference between the two densities, AA e , and their 
average value, N ave , are listed in Table U along with the 
subband energy difference A. A look at this table reveals 
the sensitivity of our analysis, which permits us to discern 
differences of less than 2 percent (sample C) in the two 
carrier populations. We can therefore say with confidence 
that in sample D, for which no double oscillations were 
observed even at the lowest fields (0.3 T), the two gases 
have essentially the same electron density. 

What the above analysis does not explain is the physi- 
cal origin of the carrier asymmetry and its difference from 
one sample to another. Could it be due to the Rashba 
effect, or is it caused by some external factor, for ex- 
ample the heterostructure growth process? Although it 
cannot be ruled out completely, it is unlikely that the 
former effect is the cause of the asymmetry summarized 
in Table I. The Rashba splitting is closely related to the 
electric field along the direction of the heterostructure's 
potential profile, which in turn is proportional to the 2D 
electron density in the accumulation potential. If the 
Rashba splitting were behind the observed asymmetry, 
then for structures of comparable carrier densities, such 
as those in Table I, one would expect similar values for 
the parameter A. In sharp contrast, A ranges from 10.6 
meV to meV. 

A far more likely explanation for the carrier asymme- 
try lies in the epitaxial growth process itself. Though the 
heterostructures were nominally symmetric with respect 
to the central GaSb layer, the actual structures might be 
slightly asymmetric. For instance, it is well known that 
dopants diffuse along the growth front. It is then possi- 
ble for the effective thickness of the undoped InAs region 
near the substrate to be smaller than that near the sur- 
face of the structure, with the corresponding asymmetry 
in the accumulation-layer profile. A small unintentional 
difference in thickness for the AlSb barriers can also re- 
sult in a difference in the amount of charge transferred to 
each of the InAs accumulation layers. Finally, the qual- 
ity of the material interfaces might affect the number of 
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interface states, which in InAs-AlSb-GaSb heterostruc- 
tures are known to be a source of extrinsic charges. 

This last possibility is vividly illustrated by comparing 
the behavior of samples C and D, which nominally have 
identical layer thicknesses. Their current-voltage charac- 
teristics were practically identical to each other from T 
= 300 K to 4.2 K, exhibiting large negative differential 
conductance even at room temperature. However, while 
sample C showed two sets of magnctoconductance oscil- 
lations [see Fig.njc)] and had A = 0.9 meV (see Table I), 
sample D showed a single set [Fig. dd)] and A = 0. In 
addition, in sample D the spin splitting was well resolved 
at a lower field. 

Samples C and D were grown one immediately after 
the other and under the same conditions, with one excep- 
tion. The growth of both samples was optimized to form 
an InSb-type bond at the AlSb-InAs interfaces, which 
has been shown to produce a smoother interface than an 
AlAs-type bondi^^iS However, sample D was prepared 
by adding an extra 1/4 In monolayer on the bottom 
AlSb-InAs interface. This extra layer reduces interface 
roughness even further, as revealed by plan-view scan- 
ning tunneling microscopy^SS 

Since the observation of a double set of oscillations 
in the magnctoconductance requires both an asymmetry 
in the carrier density and a violation of the An = 
condition, in principle it could be argued that the effect of 
the added smoothness at the interfaces of sample D is not 
to eliminate that asymmetry but just to reduce scattering 
processes that destroy Landau-level-index conservation. 
In fact, the reduction of scattering in sample D is evident 
in the appearance of ShH oscillations and spin splitting 
states at lower fields than for sample C. 

To discern between these possibilities we studied the 
tunneling magnetoconductance while a small bias was ap- 
plied between the electrodes. Such bias is not expected 



to change significantly the scattering processes but it can 
affect the potential asymmetry, as we have seen above for 
sample B. Under a 5-meV bias, the magnetoconductance 
of sample D showed a clear double-oscillation behavior 
at low fields, indicative of two different electron popu- 
lations. We can therefore conclude that the exceptional 
smoothness achieved in sample D leads to identical car- 
rier densities in the two 2D electron gases. 

In summary, we have shown that the presence of beat- 
ing and double oscillations in the tunneling magnetocon- 
ductance between 2D electron gases in nominally sym- 
metric heterostructures is due to the difference in the 
carrier densities of the two gases. The origin of this dif- 
ference lies in an extrinsic asymmetry that can be elimi- 
nated by reducing interface roughness during the epitax- 
ial growth process. In the samples with the smoothest 
interfaces the conductance showed a single set of oscilla- 
tions from the lowest magnetic fields, spin splitting was 
resolved above 2 T, and beyond 10 T pronounced fea- 
tures were observed corresponding to fractional occupa- 
tion of the magnetic levelsiSi On the one hand, these 
results cast some doubts on the apparent strength of the 
Rashba splitting and its suitability for spin devices based 
on tunneling. On the other, the results are encouraging in 
that they show that, under optimum growth conditions, 
balanced populations of high-mobility 2D gases can be 
formed in InAs accumulation layers, with which 2D-2D 
interaction effects can be probed. 
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